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Ultrasonic attenuation measurements provide nondestructive means for 
characterizing material properties. In polycrystalline solids the frequency 
dependence of attenuation characterizes microstructure and hence the mechani- 
cal properties governed by microstructure. Microstructure, morphology, speci- 
fic mechanical properties, and ultrasonic characteristics of materials are 
often uniquely interrelated (Green, 1973a, 1981, 1985; Vary, 1984). 

Although expressions that relate ultrasonic attenuation to specific 
microstructural factors exist, they are fragmented and unsatisfactory because 
they consist of diverse piecemeal descriptions and ad hoc combinations of 
attenuation mechanisms over various frequency ranges of interest (Merkulov, 
1956; Papadakis, 1960, 1965b; Serabian and Williams, 1978). A unified 
approach is needed to overcome these shortcomings and to help resolve ques- 
tions that arise in attempts to identify and evaluate key microstructural 
factors that govern mechanical properties. 

The approach given herein depends on treating material microstructures as 
elastomechanical "filters" that have analytically definable transfer func- 
tions. These transfer functions can be defined in terms of the frequency 
dependence of the ultrasonic attenuation coefficient (Vary, 1980b, 1984). The 
transfer function concept provides a basis for synthesizing expressions that 
characterize polycrystalline materials relative to microstructural factors 
such as mean grain size, grain-size distribution functions, and grain boundary 
energy transmission. Although our approach is nonrigorous, it leads to a 
rational basis for combining the previously mentioned diverse and fragmented 
equations for ultrasonic attenuation coefficients. 


BACKGROUND AND THEORY 
General 

Ultrasonic materials characterization involves both velocity and attenu- 
ation measurements. Velocity measurements are useful for measuring elastic 
constants, residual stresses, density, etc. This paper concentrates on atten- 
uation measurements because they are pivotal to understanding interrelations 
between microstructure and mechanical properties. The characterization of 
mechanical properties depends on precise measurements of attenuation. 

Ultrasonic frequency-spectrum analysis provides a means for character- 
izing microstructural factors that govern mechanical properties (Papadakis, 
1976; Fitting and Adler, 1981; Vary, 1984). The frequency range emphasized 
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herein is that normally used for interrogating polycrystalline engineering 
materials (i.e., about 10 to 100 MHz). By using ultrasonic frequency- spectrum 
analysis and signal deconvolution algorithms (Bracewell, 1978; Bendat and 
Piersol, 1980; Vary, 1980b; Williams, et al., 1985a), it is possible to define 
transfer functions that characterize microstructures relative to their effect 
on wave propagation and also on mechanical properties. It will be seen that 
these transfer functions are definable in terms of the frequency dependence of 
the attenuation coefficient. 


Pulse Echo Method 

For the purposes of this paper we assume a material sample with flat, 
parallel opposing surfaces. The sample is examined with a broadband, pulse- 
echo system using digital signal acquisition and processing (fig. 1). A buf- 
fered, piezoelectric transducer collects a set of echoes returned by the back 
surface of the sample. Only the first two echoes, Bq and B 2 , are selected 
for analysis. Signal acquisition, processing, and analysis are done as des- 
cribed by Papadakis (1976), Vary (1979b, 1980a, b), Generazio (1984), and Kline 
(1984). The echoes are transformed to the frequency domain by digital Fourier 
transform algorithms. Ultrasonic attenuation as a function of frequency and 
therefore information on material microstructure is obtained by deconvolution 
of echoes Bq and B2. 


Material Transfer Function 

In the following discussion the quantities Bq, B2, Iq, l2« T, and R 
(fig. 2) are Fourier transforms of corresponding time-domain quantities. The 
quantities Bq, B2 , Iq, and I2 are spectra of waveforms, and T is the 
material transfer function, to be defined in the frequency domain. (Phase 
spectra and such associated topics as waveform synthesis and phase interfer- 
ence effects are not considered herein.) The reflection coefficient R is 
also a function of frequency, albeit a weak function in most cases (Papadakis, 
1968; Generazio, 1984). The transfer function will be defined in terms of 
Bq, B 2 , and R. 

Internal echo Iq is the source of the "reduced" signal Bq (i.e., 
reduced by the factor (1 +■ R) , where R is the reflection coefficient at the 
specimen- transducer interface). The reflection coefficient is unity (i.e., 1) 
at the free surface. Algebraic signs carried by reflection coefficients 
depend on whether they are defined as power or amplitude quantities and on the 
input- output direction of signals across interfaces. The "power" (energy) 
definition as given in appendix A applies to R herein. A part of the energy 
of Iq is reflected and appears as the second internal echo I 2 giving the 
reduced echo B 2 ; thus 

B x = (1 - R) 1^ (1) 

B 2 = TR(1 - R)I (2) 

where Bq, B2, and Iq are the "power" (i.e., magnitude- squared spectra of 
corresponding amplitude squared time- domain signals). Combining equations (1) 
and (2) yields 
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There are transfer functions associated with the instrumentation, signal 
transduction, and other aspects of the signal acquisition system. These can 
be ignored because in the pulse-echo configuration they vanish, just as the 
term (1 + R)li vanishes in equation (3). 


Attenuation Coefficient 


Papadakis (1968; 1976) has demonstrated that the attenuation coefficient 
can be found by frequency spectrum analysis and that it may be written as 



where x is the sample thickness. The attenuation coefficient a is a func- 
tion of frequency and hence as with the quantities T, , and B 2 , it is 
appropriate to treat it as a Fourier transform of a time-domain quantity. Then 
because the transfer function T and the attenuation coefficient a. are 
both defined in terms of B^, B 2 , and R, equations (3) and (4) can be com- 
bined to get 


T = exp (- 2xot ) (5) 

Thus the transfer function is defined in terms of an attenuation coefficient. 
Conversely, if the material transfer function T can be prescribed apriori, 
appropriate expressions can be found for the attenuation coefficient ot. It 
will be seen that a general, unified expression can be derived for a. for the 
various energy loss mechanisms in polycrystalline solids by using the preced- 
ing equations and the transfer function concept. 


Attenuation Mechanisms 

The energy loss mechanisms that govern ultrasonic attenuation in the fre- 
quency ranges (bandwidths) of interest for polycrystalline engineering solids 
are known (Mason and McSkimmin, 1947; Lifsitz and Parkomovskii, 1950; Granato 
and Liicke, 1956; Merkulov, 1957; Truell, et al., 1969; Bahtia, 1976). The two 
principal causes of attenuation are scattering and absorption. Diffusion, 
Rayleigh, and stochastic (or phase) scattering losses are extrinsic; absorp- 
tion losses from dislocation damping, anelastic hysteresis, relaxation, and 
thermoelastic effects are intrinsic to individual grains (i.e., crystallites). 

We are not concerned with geometrical losses, such as those due to dif- 
fraction effects, which are not inherent to material microstructures. These 
effects can be either controlled or eliminated from attenuation measurements 
by well-known experimental and data reduction procedures (Papadakis, 1975) in 
order to get the true attenuation coefficient as a function of frequency. 
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Extrinsic Mechanisms 


Scattering usually accounts for the greatest portion of losses in poly- 
crystalline aggregates (Papadakis, 1964a, b, 1965a, b, 1970; Serabian and 
Williams, 1976; Goebbels, 1980). For polycrystalline aggregates expressions 
that have been developed for the scatter attenuation coefficient a as a 
function of frequency f are referred to three regimes that are defined by 
the ratio of mean grain size a to dominant wavelength X. 

For the Rayleigh scattering regime, where X >> ira, 


3.4 

a = c a f 
r r 


For the stochastic (phase) scattering regime, where 


,2 

a = c af 
P P 

For the diffusion scattering regime, where X « ira, 


-1 

a d = c d a 


X s ira, 


( 6 ) 


(7) 


( 8 ) 


In equation (6) a^ represents grain volume, where the grains are to be con- 
sidered as having a variety of nonequiaxed, nonspherical shkpes (Papadakis, 
1965b). The constants c contain other quantities (i.e., geometric factors, 
longitudinal and transverse velocities, density, and elastic anistropy fac- 
tors) (Mason and McSkimmin, 1947, 1948; Lifsitz and Parkomovskii, 1950; 
Merkulov, 1957; Truell, et al., 1969). These factors will be considered later 
in this paper. 


Intrinsic Mechanisms 


Absorption losses due to dislocation damping, hysteresis, thermoelastic 
effects, etc., are intrinsic to grains (crystallites) and involve direct con- 
version of acoustic energy to heat. These attenuation mechanisms are essen- 
tially independent of grain size, shape, volume, etc. 

For hysteresis losses (Kolsky, 1953; Nowick and Berry, 1972) 


Ot = c f 
h h 


For thermoelastic losses (Lucke, 1956) 


- °t f 


(9) 


( 10 ) 


Hysteresis losses (eq. (9)) arise when acoustic waves cause stress-strain 
damping. Hysteresis losses with a first-power frequency dependence are usu- 
ally observed in single crystals and amorphous solids and with difficulty in 
polycrystalline solids (Mason, 1958). Equation (10) gives the frequency 
dependence of thermodynamic losses that arise when longitudinal waves produce 
heat flow from dilatation to compression regions. 
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Viscous losses (Auld, 1973) also exhibit a second-power frequency depen- 
dence but are generally negligible in solids. Models for absorption losses 
caused by dislocation vibrations, relaxation effects, and internal friction 
predict second-power frequency dependence down to frequency independence 
(Granato and Lucke, 1956; Seeger and Schiller, 1962; Truell, et al., 1969). 
Absorption losses due to magnetic domains, electrons, or phonons are special 
cases involving ferromagnetic materials, very high frequencies, or cryogenic 
temperatures, respectively, and are not considered here. 


Combined Expressions 

Total attenuation coefficients are usually written ad hoc as sums of 
coefficients for scattering and absorption. For example, equations (7) 
and (9) are combined to form an expression for the total attenuation coeffi- 
cient a, where 


% - V * y f2 (11) 

or equations (6) and (9) are combined to get 

3 4 

a = c. f + c a f (12) 

hr h r 

Combined expressions as in equations (11) and (12) are convenient for fitting 
experimental data and for analyzing the contributions of attenuation mecha- 
nisms and underlying microstructural factors. 


An inconsistency occurs in attempting to combine the expressions for 
Rayleigh and stochastic attenuation in equations (6) and (7) because as ultra- 
sonic frequency increases, the fourth-power term for Rayleigh scattering 
becomes dominant. This contradicts the expectation that for higher frequen- 
cies, as wavelength diminishes and approaches grain size, the second power of 
frequency due to stochastic scattering should dominate. Papadakis (1965b) and 
others have dealt with this issue by applying equations (6) and (7) piecemeal 
and adapting a smoothing function to form continuity between the Rayleigh and 
stochastic regimes. Smith (1984) has observed that the transition from 
Rayleigh to stochastic scattering depends on grain-size distribution and sug- 
gests introducing a probability function in scatter attenuation equations pro- 
posed by Roney (1950) and Serabian and Williams (1978). 

Maisel (1984) overcomes the problem of expressing the transition from 
Rayleigh to stochastic scattering by using an electrical network analog to 
model ultrasonic scattering. Maisel’ s results suggest expressions for the 
attenuation coefficient that combine Rayleigh and stochastic frequency depen- 
dencies in a continuous function. One such expression is 


a 

rp 


c 

rp 


4aV 

1 + 4a 2 k 2 


(13) 


where c is a material constant incorporating average fractional deviations 
of density and compressional modulus, a is mean scatterer size, k is wave 
number (2irf/v), and v is velocity. The formulation of equation (13) allows 
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the fourth power of frequency to dominate at low frequencies and the second 
power to dominate at high frequencies. A similar formulation was derived by 
Perkeris (1947) by assuming an autocorrelation function for polycrystalline 
inhomogeneity. Although both formulations allow a smooth transition from 
fourth- to second- power frequency dependence with increasing ultrasonic fre- 
quencies, no allowance is made for frequency independence at extreme frequen- 
cies in accordance with equation (8). Rokhlin (1972) describes a generalized 
formulation that encompasses Rayleigh, stochastic, and diffusion scattering. 
However, in Rokhlin* s approach a smooth transition between scattering regimes 
is lost because of the periodic nature of the sine function that was used. 


MICROSTRUCTURE CONSIDERATIONS 
Grain Size 

Before dealing with transfer function synthesis we briefly review some 
relevant aspects of the interplay of microstructure and the previously dis- 
cussed formulations for ultrasonic attenuation. 

A universally cited quantity for characterizing polycrystalline micro- 
structures is "mean grain size.** This quantity is used regardless of diffi- 
culties inherent in measuring it or assigning appropriate values to it, 
especially in materials that exhibit complex microstructures (subgrains, 
second phases, etc.). Nevertheless it is appropriate, for many single-phase 
polycrystallines with "uniform" microstructures, to speak of a mean grain 
size. This depends on discovering a definitive grain-size distribution func- 
tion (DeHoff, 1968). Most polycrystallines exhibit grain-size distributions 
that are highly skewed and nonGaussian because of a high incidence of large 
grains (Papadakis, 1964b, 1965b; Schuckler, 1968; Evans, et al. 1978; Smith, 
1984; Tittmann, et al., 1985). 

There are polycrystalline aggregates that have microstructures with 
readily defined mean grain sizes (Schuckler, 1968). When expressions of the 
form of equations (11) and (12) are fitted to experimental attenuation data 
for samples of these materials, a new expression for the same sample after 
annealing to a different mean grain size can be obtained by modifying the 
coefficients of frequency by a numerical scale factor (e.g., the ratio of mean 
grain sizes). This has been demonstrated for polycrystalline samples of 
nickel and copper where measurements were predominantly in the stochastic 
scattering regime (i.e., principal wavelengths commensurate with grain size 
(Generazio, 1985)), A general scaling parameter for grain size has been 
described by Serabian and Williams (1978). 

Scaling of attenuation expressions probably applies only to homomorphic 
microstructures having grain-size distribution functions that are Invariant 
with changes in mean grain size. Photomicrographs of the successive heat 
treatments produced by Generazio (1985) seem to scale (i.e. , map) exactly in 
this sense. Apparently scaling can be done even with rather complex micro- 
structures despite difficulties in producing or identifying homomorphic grain- 
boundary patterns. 
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Scattering Centers 


Scattering theories usually consider ensembles of scattering centers 
embedded in a featureless continuum. Current theories for polycrystallines 
account only for scattering by equiaxed grains and neglect effects of texture, 
anisotropy, grain-size variations, etc. Indeed, well-formulated theories for 
attenuation exist only for the simple case of single-phase polycrystalline 
solids with virtually identical, equiaxed grains and for frequencies that 
satisfy the conditions ka << 1 or ka >> 1, where k is wave number and a 
is mean grain (scatterer) size. 

Underlying equations (6) and (7) is the assumption of single scattering 
by independent scatters; the conceptual and theoretical difficulties that 
arise with multiple scattering where scatterers interact and velocity disper- 
sion effects are likely to arise are ignored. For example, equation (6) is 
based on the "Rayleigh criterion" and ideal elastic scatterers that are suffi- 
ciently separated to ensure single scattering (Truell, et al., 1969). 

Fourth- and second-power relations given in equations (6) and (7), 
respectively, have been experimentally confirmed for special cases (Winkler 
and Murphy, 1984). Papadakis (1981) and Sayers (1985) realized the fourth- 
power relation of equation (6) by measuring scattering due only to "minority" 
phases (e.g., carbon nodules and inclusions) in polycrystalline aggregates. 
Evans, et al., (1978) used a similar approach for materials with sparse dis- 
tributions of "predominant" scatterers in a polycrystalline "continuum." 

Apparently, when used with materials exhibiting discrete, widely sepa- 
rated, "strong" scattering centers, the Rayleigh criterion is satisfied and 
the difficulties due to multiple scattering are avoided. It seems remarkable 
that frequency dependences predicted by rather simplified scattering models 
can be realized even with complex polycrystallines . 


Empirical Correlations 

For most polycrystallines the overall experimentally determined frequency 
dependence of attenuation is not an integral power. Exponents on frequency 
even in the Rayleigh or stochastic regimes rarely equal exactly 4 or 2, 
respectively (Papadakis, 1981). This is probably due partly to imprecisions 
in experimental measurements and to attenuation data being defined only within 
the bandwidth used to interrogate the sample. Within limited bandwidths the 
best fit of experimental data simply reflects only a portion of the actual 
attenuation function. 

Because of the aforementioned factors it is often possible to fit attenu- 
ation data with expressions that allow the exponent on frequency to be an 
experimentally determined variable (Vary, 1980a, b). Among the simplest 
expressions of this kind is 


>a< = cf m f x < f < f 2 (14) 

where denotes that a is empirically defined only in a limited fre- 

quency range and where m may assume noninteger values. The quantities c 
and m will change correspondingly with microstructural changes that affect 


263 



attenuation (mean grain size, grain size distribution, etc.). Experimental 
results show that equation (14) is appropriate for fitting attenuation data 
over a frequency range that spans the Rayleigh scattering regime for a variety 
of polycrystalline materials (Papadakis, 1964a, 1976; Vary, 1978; Bozorg- 
Grayeli, 1981; Vary and Hull, 1982, 1983). It will be seen that equation (14) 
gives statistically valid descriptions of data in the Rayleigh regime. 


Grain Interfaces 

From a purely physical standpoint grain facets and surface areas should 
influence mechanical properties of polycrystallines. This is certainly true 
of properties that depend on the surface energy of grains, properties affected 
by the grain-boundary thickness, properties for which grain boundaries are 
obstacles, and properties connected with grain- and phase-boundary migrations 
and obliterations. The transfer function concept, described next, provides a 
basis for a unified expression for attenuation that incorporates the previ- 
ously discussed attenuation mechanisms and also a means for studying the 
effects of variations in grain-boundary energy transmission factors. 


TRANSFER FUNCTION SYNTHESIS 
Approach 

An attenuation transfer function for polycrystalline solids is defined in 
terms of a lattice model that incorporates the effects of grain-boundary 
(interface) spacing and energy transmission across interfaces. The lattice 
model is generalized by the introduction of probability distribution functions 
to account for variations in grain and scatterer size and boundary spacing. 

The transfer function concept provides a convenient and rational basis for 
incorporating the previously described attenuation mechanisms into a unified 
expression for ultrasonic attenuation. 


Interface Transfer Function 

The transfer function of an interface between two media is defined in 
terms of the reflection coefficient at the interface. We assume ideal flat 
interfaces between dispersion-free media with unequal acoustic impedances. In 
this case the reflection coefficient R is taken as a frequency- independent 
Heaviside step function in the frequency domain, corresponding to a delta 
(impulse) function in the time domain, see figure 3 and appendix B. 

The reflection coefficient R in equations (1) to (4) conforms with the 
interface transfer function notion given here. If it becomes necessary to 
consider velocity dispersion or nonideal (thick, strained, acoustically rough, 
etc.) interfaces, expressions for R (and the quantity G, defined below) 
would simply assume frequency-dependent forms (Papadakis, 1976; Kinsler, 
et al., 1982; Generazio, 1984). 
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Lattice Model Transfer Function 


The lattice model illustrated in figure 4 represents a volume that con- 
sists of series of plane parallel interfaces separating thin layers with 
unequal acoustic impedances (compare with Rokhlin (1972), where a volume is 
divided into equal rectangular parallepipeds) . The model is one dimensional 
and accounts for energy losses suffered by longitudinal rays traveling at 
normal incidence to the interfaces. Each of the n internal interfaces has a 
transfer function defined in terms of an internal reflection coefficient G 
(appendix B) . For the moment consider each layer as a featureless (scatterer 
free) medium with an intrinsic attenuation function 

A = | A( f ) | = exp (~a^8,) (15) 

where 8. is layer thickness and A is the magnitude (power) spectrum of 
A(f ) . As developed in appendix C the transfer function T for the lattice 
model can be reduced to the form 

T = (H - G) 2n A 2(n+1) []C. . .] (16) 

where H = |H(f)| =1 is the Heaviside unit step function that will be ren- 
dered as 1 hereinafter (i.e., (H - G) = (1 - G)). The multiple convolution 
(multiple reflection) factor [£...] in equation (16) is 



n n 

2 

n - 

r 1 

rE...] = 

£ (RG)' 1 -V (<1 - 1) 


]T G q_1 A 2(q_1) 

£ A'VH-U 


_1 


_i 

.1 


2(n-l) 


(17) 


The number of first-order multiple reflections in each lattice layer is 
q - 1, and the number of layers is n + 1. Equation (16) defines T in terms 
of power (energy) reflection functions R and G (appendix A). The terms 
A, G, R, and T are Fourier transforms. The various powers and products of 
A, G, and R are convolutions of their time-domain counterparts. For exam- 
ple, RG in equation (17) is the convolution of corresponding time-domain 
quantities. The quantity []^...] given in equation (17) explicitly contains 
only first-order multiple reflections. It will be seen that this is suffi- 
cient for generalizing the model. 


Lattice Model with Scattering 

The lattice model may be extended to include scattering by taking A in 
equations (16) and (17) to consist of both intrinsic and extrinsic attenuation 
factors , 


A = [exp(-a. 8.) ] [exp(-a 8.)] (18) 

l e 

If it is postulated that each layer contains an independent scattering center 
in addition to intrinsic attenuation, equation (16) takes the form 

T = (1 - G) 2n [exp(-a i 8.)exp(-a e SO] 2(n+1) [i;...] (19) 
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Let the intrinsic attenuation be due to anelastic hysteresis, and in accordance 
with Nowick and Berry (1972) express it as 


a. 

l 


" C h f 



( 20 ) 


where h is a hysteresis constant and v is velocity. For extrinsic attenu- 
ation take the Rayleigh scattering expression derived by Mason and McSkimmin 
(1947), 


a 
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3 C 4 
c a f 
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(l6v\\ 3„4 

llT/ af 


( 21 ) 


where a is mean scatterer size (with 4ira /3 replacing the mean grain vol- 
ume in the Mason and McSkimmin (1947) derivation) and K is elastic anisot- 
ropy (i.e., average change in elasticity encountered from one layer to the next 
(Mason and McSkimmin, 1948; Lifsitz and Parkomovskii, 1950; Merkulov, 1957)). 


Generalization for Polycrystallines 

We generalize the lattice model to resemble a polycrystalline by postu- 
lating unequal layer thicknesses and nonparallel interfaces. Thus any ray 
traveling at normal incidence to the front and back surfaces (still mutually 
parallel) encounters a sequence of interfaces similar to those in a polycrys- 
talline aggregate. 

Combining equations (4) and (19), noting that, for large values of n, 
(n/x) ■- (1/8.) while (n + 1) 8, = x and taking 8. as mean layer thickness, an 
expression is obtained for a in terms of T, 

a = -Qj In (1 - G) + + a e + In [£...) (22) 

The quantity G is now taken as average reflection coefficient for the ran- 
domly oriented interfaces. If G is sufficiently small and energy is indeed 
scattered away in lattice layers, the contribution of higher order terms in 
the summations in equation (17) may be ignored. This is equivalent to setting 
q = 1 in equation (17) so that In [£...] = 0 (see DISCUSSION). Although 
the energy transmission coefficient 1 - G may approach unity, the interface 
transparency term In (1 - G) is always finite (i.e., a nonzero negative). 

Given the preceding argument and noting that ft, corresponds to a mean 
grain size represented by the quantity a and combining equations (20) to 
( 22 ), 

a = -(j) In (1 - G) + (j)f + asV (23) 

Although a spherical scatterer with volume 4ira^/3 is implicit to equa- 
tion (21), in equation (23) scatterer geometry is made abstract so that s is 
an effective mean scatterer diameter while a is mean scatterer length. Thus 
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we allow for nonequiaxed grains or scatterers with aspect ratios a/s of 1 or 
greater. Rokhlin (1972) assumes grains to be parallelepipeds having square 
ends but different side lengths in the direction of wave propagation. 


Generalization for Grain-Size Distribution 

The range of grain sizes usually contained in a polycrystalline can be 
described by a histogram or, for a large grain population, by an analytical 
probability distribution function. Grain-size distributions are usually 
nonGaussian and their probability distributions are typically log-normal func- 
tions (DeHoff, 1968; Schtickler, 1968; Smith, 1984). A function that describes 
grain-size distributions representative of most polycrystalline solids is 
given by 


p(y) 



p(y) = 0 for y < 0 


(24) 


where y is a spatial (length) variable and a is mean grain size. We have 
pointed out that the quantity a is defined in the frequency domain and can 
therefore be treated as the Fourier transform of a time-domain quantity. To 
incorporate probability distribution functions in the expression for a given 
by equation (23), we transpose to the time domain, take the Fourier transform, 
and find the results of convolution with probabalistic factors in equa- 
tion (23). Therefore let 



p(t) = 0 for y < 0 


(25) 


where u is mean transit (grain delay) time corresponding to mean grain size 
a and a = uv, where v is average velocity (fig. 5). The Fourier transform 
of p(t) is 


P(f) = jjL + (iruf)^J exp£~2i tan ^(iruf)J (26) 

The phase spectrum given by the exponential term in equation (26) is ignored 
since we are interested only in the magnitude spectrum of P(f) (i.e., 

P = | P ( f ) | ) . The implicit assumption is that phase interference effects are 
insignificant. In general, the probability distribution function of the scat- 
terer diameter s is not identical to that of its axial length a. Therefore 
we admit two separate expressions for P, 


p a ' t 1 + <"V> 2 ] 
p s ‘ [ l * <*V> 2 ] 


(27) 

(28) 


where subscripts s and a refer to mean scatterer diameter and length, 
respectively. That is, u a = a/v and u s = s/v. 
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Now consider the first term in equation (23) as a probabilistic quantity. 
In the frequency domain the probabilistic weighting of the reflection function 
for transit time variations is GP and the probabilistic weighting for trans- 
mission is G(1 - P) . Therefore the frequency dependence of the boundary 
energy transmission term 1 - G in equation (23) is In [1 - G(1 - P)]. The 
expression for P in equation (27) will be used for this case. The middle 
term in equation (23) is independent of grain size and remains unaltered. In 
the third term of equation (23) a and s are probabilistic so that 


a 



In [1 - G(1 - P )] + 
a 



/ l6ir 4 K \ 
\ 3vU a ) 


2 2 4 _ 

u u f P P 
as as 


(29) 


where u a = a/v, u s = s/v, and v is average velocity. The corresponding 
probability functions P a and P s are to be defined according to the nature 
of the various microstructural factors that govern attenuation (grains, pre- 
cipitates, etc.) in a particular solid. A further generalization would take 
grain-boundary spacing and scatterer spacing (and dimensions) to be indepen- 
dent of each other. 


VERIFICATION AND APPLICATIONS OF MODEL 


Comparison with Scattering Theory 


We examine a special case of equation (29) wherein grains are equiaxed 
scatterers (i.e., a = s = vu) . Also, take P a = P s = P = f P ( f ) | . For the 
moment ignore the first two terms in equation (29) so that consequences of the 
preceding derivations may be compared with classic expressions for scattering. 
The extrinsic attenuation (third) term of equation (29) expands to 
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e 


16it K 

V 3vu ) 


iml 

[l + (iruf) 2 ] 


(30) 


For \ » ira equation (30) reduces to the expression for Rayleigh scattering 
(Mason and McSkimmin, 1947), 
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(31) 


For \ = ira equation (30) reduces to the expression for stochastic (phase) 
scattering, 
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ep 




af 


(32) 


For \ << ira equation (30) reduces to the expression for diffusion scattering. 


a 


ed 



(33) 
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Thus we have in equation (30) an expression for extrinsic attenuation that 
covers all theoretically postulated scattering mechanisms continuously over 
all frequencies. 

Although we have adopted the scatter attenuation expression of Mason and 
McSkimmin (1947), the form developed by Lifsitz and Parkomovskii (1950) would 
be appropriate for defining extrinsic attenuation in equations (21) and (30) 
to account for elastic randomization and attendant effects of mode conversion. 
However, it is convenient and sufficient for our purposes to continue using 
the expression for extrinsic attenuation given in equation (21). 


Comparison with Experiment 


Consider a simple polycrystalline where, as in the previous section, 
grains and other scatterers are taken to be identical, so that a = s = vu, 
where a is mean grain size, u is mean grain delay, and v is mean veloc- 
ity. Using equation (30) for extrinsic scattering, equation (29) becomes 
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In 
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1 + (iruf) 
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(34) 


We assign appropriate values to G, K, and h for comparing attenuation 
coefficients predicted by equation (34) with experimental data for nickel, 
copper, and iron polycrystallines. 


Serabian and Williams (1978) estimated hysteresis constants for several 
metals (e.g., h = 2x10” 5 f or 97,4 percent iron alloys and h = 3xl0” 4 for 
98.9 percent nickel alloys). Herein we assume h = 3x10”^ for 99.99 percent 
copper. 


Both G and K are estimated from equivalent equations and may be taken 
to be identical for cubic crystals (Mason and McSkimmin, 1947, 1948; Mason, 
1958): 



where c' is the elastic tensor coefficient (compression modulus) and <c*> is 
its average value. The z = pv terms are defined in appendix A and are 
based on principal longitudinal velocities. We can assume that for any given 
material G is directly proportional to K and that for cubic crystallites 
G = K. According to Mason and McSkimmin (1948) G = K « 0.0074 for copper, 
0.0067 for iron, and 0.0072 for nickel. 

Attenuation-versus-frequency curves generated with equation (34) are 
compared with experimental curves for the cubic materials iron, copper, and 
nickel in figures 6 to 8. The estimated values for K and G do not account 
for the effects of oblique incidence or mode conversion at boundaries; these 
effects would increase their effective values. Therefore the values for K 
cited above should give low predictions of attenuation since shear waves pro- 
duced by mode conversions are more effectively scattered. Elastic anisotropy 
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for shear waves can exceed that for longitudinal waves by perhaps a factor of 
10. It is probably for this reason that the best agreement with experimental 
data is obtained when we double the estimated values for K for use in equa- 
tion (34). It will be seen that K is in general a composite quantity 
containing other scattering factors. 

Figure 6 compares original experimental data of Merkulov (1957) for iron 
with attenuation-coefficient-versus- frequency curves generated with equa- 
tion (34). Excellent agreement between the experimental and predicted attenu- 
ation coefficients is evident in figures 6(a) and (b) for the particular values 
selected for K, G, and h for iron, even to the extent that curve crossovers 
at approximately 14 MHz are predicted. Figures 6(c) and (d) exhibit attenu- 
ation curves for iron over a broader frequency range of 100 MHz, as opposed to 
the 20- MHz range in figure 6(b). Figure 6(c) is based on K = 0.0067 and 
figure 6(d) on K = 0.0134 in equation (34) for grain sizes from 40 to 250 ym. 
Figures 7 and 8 compare original experimental data of Generazio (1985) for 
copper and nickel, respectively, with attenuation- versus- frequency curves 
generated with equation (34). 


Data Simulation for Curve Fitting 

It is apparent that prior knowledge of the material parameters h, G, 
and K is needed to predict the attenuation coefficient with equation (34). 
Estimates of these quantities were sufficient to simulate attenuation curves 
for iron, copper, and nickel in figures 6 to 8. We illustrate an empirical- 
graphical approach for estimating and verifying values for h, G, and K in 
solids, where these quantities are not easy to determine. 

Our approach is founded on the observation that over a limited frequency 
range attenuation data can often be fitted quite well with equation (14) 
(Papadakis, 1964a; Vary, 1978, 1982; Bozorg- Grayeli , 1981). It has been shown 
that for data in the Rayleigh regime extrapolations based on the empirical 
constants c and m lead to excellent correlations with mechanical proper- 
ties (e.g., fracture toughness) (Vary , 1978; 1979; Vary and Hull, 1982, 1983). 

Ho analytical proof is offered to show why an equation of the form of 
equation (14) should be expected to represent either equation (34) or (29) in 
any particular frequency regime. However, we will show that equation (14) 
satisfactorily represents experimental data in restricted frequency regimes by 
using equation (34) to generate raw data sets in the Rayleigh regime. The 
point is illustrated by simulating data sets for a range of mean grain size 
a and values of K, G, h, and v that represent iron. 

A series of discrete data sets generated with equation (34) are exhibited 
in figure 9 for the frequency range 10 to 100 MHz, where the Rayleigh wave- 
length criterion is satisfied (i.e., X >> ira) . We ignored data below 15 MHz, 
where diffraction effects are most pronounced for actual data sets. Treating 
the simulated data as actual raw attenuation data, we determined c and m 
values by a least- squares fit to equation (14) for each “material sample.*' 

The c and m values and the correlation coefficients of the curve fittings 
were tabulated for each data set. As seen in figure 9 calculated frequency 
exponents m were clearly nonintegral and ranged from approximately 1.8 to 
2.3, corresponding to mean grain sizes from 2 to 20 ym, in the case illus- 
trated. The simulated exponents spanned a range similar to experimentally 
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determined frequency exponents for a number of polycrystallines in the Rayleigh 
regime (Vary, 1978). 

For the material parameters chosen for generating data for figure 9, the 
correlation coefficient was always greater than 0.999. Thus by treating the 
simulated data precisely the way one could treat a set of actual attenuation 
data, we found that equation (14) renders a statistically valid representation 
within the frequency bounds specified. 


Interrelation of Parameters 

In addition to representing attenuation data for a given material sample 
within specified frequency bounds, the parameters c and m are interrelated 
for variously conditioned samples of the same material. As an example, exper- 
imentally determined values of c and m for two maraging steels were plotted 
and tabulated (fig. 10), All materials for which the quantities c and m 
have been experimentally measured exhibit the strong functional relation 
between c and m evident in figure 10. The quantities c and m are evi- 
dently interdependent for any set of related materials or samples of the same 
material that apparently differ due to metamorphic changes brought about by 
heat treatment or other thermomechanical processing methods that retain global 
microstructural patterns but alter mechanical properties. 

By varying the material parameters G, K, h, and v in equation (34) and 
allowing the mean grain size to take representative values, the c and m 
quantities corresponding to any particular material can be duplicated. Para- 
metric fields based on varying G, K, and h appear in figure 11. The fields 
in figure 11 were created by using the curve- fitting method illustrated in 
figure 9 (i.e., by simulating attenuation data with equation (34), fitting 
with equation (14) in the 15- to 100- MHz range and including in the fields 
only points for which the curve fit correlation exceeded 0.999). Values 
selected for h, G, and K corresponded to those for maraging steel and v 
was taken to equal 0.56 cm/ys, the mean value for a 250-grade maraging steel 
(Vary, 1978). 

Coplotted in each field in figure 11 are data for the 250-grade maraging 
steel from figure 10 for which the mean grain size ranged from about 8.5 to 
13 ym. Accordingly, in figure 11(a) the maraging steel data fall near the 
curve for mean grain size a equal to 10 ym. This is not fortuitous, since 
the location and shape of the field, relative to experimental c and m 
values for the maraging steel, depend on judicious selection of values for h 
and G. Figures 11(b) and (c), where we assume that a = 10 ym, appear to 
confirm the values G = 0.007 and h = 0.00002 for the 250-grade maraging 
steel . 

An unexpected result is that each field in figure 11 indicates that the 
quantity K, the elastic anisotropy factor, varies by approximately a half 
order of magnitude among the maraging steel samples represented. Similar 
results were obtained for the 200-grade maraging steel data in figure 10. 

This supports the notion that a material will have a different value of K 
after thermomechanical conditioning. Indeed, it has been suggested that K 
should be viewed as a composite quantity that includes a scattering coefficient 
(Mason, 1958; Papadakis, 1965a; Serabian and Williams, 1978). Heat treatment, 
strain, deformation, etc., alter attenuation properties due to effects on 
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dislocations and grain boundaries (Serabian, 1985). These effects are apparent 
when experimental values for c and m are coplotted in appropriate para- 
metric fields as in figure 11. 


DISCUSSION 

General 

Our goal has been to develop a unified, comprehensive expression for 
ultrasonic attenuation. The incentive was to find a rational basis for iden- 
tifying and evaluating key microstructural factors that govern mechanical 
properties. We noticed that although previous expressions for attenuation are 
theoretically self-consistent, they provide only piecemeal, fragmented des- 
criptions of ultrasonic attenuation mechanisms over the broad frequency range 
relevant to materials characterization. 

We describe a transfer function concept and a heuristic lattice model 
that form bases for incorporating well-known attenuation mechanisms into a 
unified expression for ultrasonic attenuation. Although we ignore phase 
interference effects, our generalized analytical expression continuously 
accounts for the waxing and waning of all of the attenuation mechanisms pre- 
viously embodied in separate equations. The generalized equation is shown to 
predict salient features of experimental attenuation-versus-f requency curves 
obtained for typical polycrystalline solids. 

We believe that our chief contribution is the transfer function approach, 
in which expressions for attenuation coefficients are recognized and handled 
as frequency- domain transforms of time-domain quantities. 


Microstructure Encompassed 

We are concerned with polycrystallines, where ultrasonic attenuation 
depends primarily on such extrinsic microstructural factors as grain- size dis- 
tribution and properties of grain interfaces. Using equation (34) we concen- 
trate on relatively simple fine-grained, isotropic, homogeneous materials with 
mean grain sizes in the range 1 ym to >100 ym. We expect that the concepts 
given herein apply equally well to anisotropic materials and to amorphous or 
single -crystal materials that exhibit microporosity, dispersoids, etc. For 
example, equation (29) admits additional variables to account for grain ani- 
sotropy and bimodal distribution functions for dispersed scatterers. 


Probability Distribution Functions 

It should be apparent that our selection of a probability distribution 
function for use in equation (29) and equations developed from it was quite 
deliberate. In addition to representing typical polycrystalline grain popula- 
tions, the probability distribution function given in equation (25) has salu- 
tary mathematical properties. These are evident in the consequences exhibited 
by equations (31) to (33), which depend on the form of the probability func- 
tion. Moreover, once the mean is specified for the distribution function in 
equation (25), the variance is specified also (i.e., the variance equals half 
the square of the mean). Alternative distribution functions, such as Gaussians 
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or combinations of normal and log-normal, may provide better descriptions of 
particular grain populations. Our immediate purpose was to indicate a rational 
method for introducing probabilistic factors in attenuation equations. 


Interface Transfer Function 

We applied the transfer function concept to a generalized version of the 
lattice model to represent polycrystalline materials. One consequence of the 
model was the introduction of a grain- boundary (interface) energy transmission 
function (e.g., represented by the term (1/a) In (1 - G) in eq. (23)). This 
term is similar to the diffusion scattering expressions of equations (8) and 
(33). It is related to K, the elastic anisotropy factor in equation (33). 
Unlike the extrinsic (scatter) attenuation terms that include K, the 
In (1 ~ G) term does not vanish in the absence of scatterers. Our model 
requires (and retains) it to account for the energy lost by reflection at 
internal interfaces. In general, G is frequency dependent (as in the case of 
nonideal, thick, or acoustically rough interfaces, see appendix A). 


Multiple Reflection Term 

In our generalization the multiple reflection term of equa- 

tion (17) is ignored. This is justifiable if the model represents polycrys- 
tallines that exhibit no significant forward- or backscatter. The waveform 
shown in figure 1 is a typical of this case (i.e., where no evidence of scat- 
ter wavelets appears either between or near the main echoes). This would not 
be true for coarse-grained materials that can exhibit considerable backscatter. 
In general, the. [£...] term needs to be considered (e.g., when first- and 
second-order multiple reflections are significant as in layered structures 
(composite laminates)). 

The quantity [£...] collects all multiple reflections that are not 
attenuated or deflected out of the interrogating beam. Referring to equations 
(17) and (22), we notice that, in our model, as G -> 0 the lattice interfaces 
become transparent, and overall attenuation will be due to intrinsic absorp- 
tion and extrinsic scattering, say, due to point scatterers, voids, etc. In 
this case equation (22) reduces to <* = + a e . Note that for all values of 

G the first term of each summation in equation (17) equals 1. 

When the average internal reflection coefficient G lies between 0 and 
1, the interfaces assume varying degrees of opacity. For sufficiently large 
values of G we can no longer ignore the first and higher power terms associ- 
ated with q = 2, 3,..., in equation (17). In this latter case backscatter 
would become significant and the effects of the reflection coefficient R at 
the transducer- specimen interface enter the expression for a, as for example 
in equation (22). Under this condition as in the case of equation (4) meas- 
urement of a depends on apriori knowledge or independent measurement of R. 


Power Functions of Frequency 

Integral exponents on frequency, as given by equations (6) to (10), should 
not be expected to apply to experimentally determined relations between the 
attenuation coefficient and frequency. Our analysis and experimental evidence 


273 



show that even over narrow frequency ranges attenuation is a nonintegral power 
function of frequency. Exponents on frequency will assume values between 1 
and 4 according to the degree of influence of material parameters such as 
those appearing in equations (29) and (34). Indeed, figures 6 to 8 show that 
exponents are themselves functions of frequency. 

Heretofore, we and other investigators have found it practical to use 
nonintegral, fixed exponents to describe empirical relations between attenu- 
ation and frequency within bounded regimes. By using the generalized attenu- 
ation equations we show that this approach is statistically valid (e.g., as in 
the case of fig. 9). However, it is now clear that considerably more informa- 
tion concerning microstructure may be extracted by applying the generalized 
expressions for attenuation in equations (29) and (34). 


Elastic Anisotropy and Related Factors 

Our preliminary analysis of the implications of the interrelation of the 
empirical quantities c and m from equation (14) confirmed the importance 
of considering variations in the effective values of K (and G) in applying 
equation (34), especially to polycrystallines modified by thermomechanical 
processing. We infer that changes can occur in the effective values of K 
that transcend changes in grain geometry (size, shape, etc.). These corre- 
spond to changes in mechanical properties (fracture toughness, yield strength, 
etc.) in materials in which only slight or no changes in grain size have 
occurred. In effect, equations (29) and (34) appear to provide entry to anal- 
ysis of subtle metamorphic changes that arise from thermomechanical processing 
of polycrystallines. 


CONCLUDING REMARKS 

A unified expression for ultrasonic attenuation has been introduced to 
supplant currently diverse and fragmented equations. A transfer function con- 
cept and a heuristic lattice model with an interface energy transmission factor 
have been described as rational bases for combining a variety of attenuation 
mechanisms and probabilistic factors in a comprehensive equation. Although our 
model is essentially one dimensional and we ignore phase interference effects, 
the resultant generalized equation predicts salient features of experimental 
attenuation-versus-f requency curves obtained for typical polycrystalline sol- 
ids. Some consequences of the unified expression have been explored relative 
to characterization of key microstructural factors that govern mechanical 
properties of polycrystallines. 
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APPENDIX A 


REFLECTION COEFFICIENT AND ACOUSTIC IMPEDANCE 


The reflection coefficient can be defined either in terms of power (energy 
intensity) or amplitude (pressure) (Krautkramer, 1969; Szilard, 1982). Taking 
R as the reflection coefficient and Q as the transmission coefficient across 
an interface, from medium 1 to 2, we have for the power case, 


and for the amplitude case, 


R + Q = 1 
P P 
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(4Z 2 Z 1 ) 
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where z^ and z ^ are acoustic impedances of media 1 and 2, 


(A3) 


p is density, and v is velocity. 

In polycrystal lines grain- boundary interfaces will reflect as well as 
transmit ultrasonic energy because of acoustic impedance "jumps" across dif- 
ferently oriented contiguous crystals. 

The preceding equations for reflection and transmission apply to ideal 
interfaces that have no thickness. If there is an intervening medium that has 
a finite thickness d such that d greatly exceeds the ultrasonic wavelengths 
transmitted, d » X, the transmission and reflection coefficients are fre- 
quency dependent (Kinsler, et al., 1982). Letting V3 be velocity in the 
intervening medium (medium 3) between media 1 and 2 and Z3 be its acoustic 
impedance, R becomes a function of frequency f with the magnitude 


z = p v 
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When ^3 = (v3/2irf) » d, equation (A4) reduces to equation (Al) 
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APPENDIX B 


REFLECTION COEFFICIENT TRANSFER FUNCTION 

In this paper the reflection coefficient R is always taken as a 
frequency-domain quantity. For an ideal interface between two nondispersive 
media with unequal acoustic impedances, R is independent of ultrasonic fre- 
quency (Papadakis, 1976; Generazio, 1984). We define R as identical to the 
Heaviside step function H(f), (Bracewell, 1978) (fig. 3) 

I R( f ) I = |H(f)| (Bl) 

where 

! 1 for f > 0 

(B2) 

0 for f < 0 

The time-domain counterpart (i.e., inverse Fourier transform) of R(f) is then 
the delta (impulse) function 6 defined by 



6 (t) dt = 1 


(B3) 


where 6=0 for t ^ 0 and t is time. 

We adopt the additional formalism that the delta function is the inverse 
transform of H(f) with arbitrary magnitude and phase spectra. This allows 
the delta function to be positioned at any point in the time domain and to 
have any value w (Bendat and Piersol, 1980) so that 



where 


for t = to 

(B5) 

0 for t * to 

and where |H(0) | = w by definition of the Fourier transform of 6(t). In 
a nondispersive medium |H(f) | = w, a constant over all frequencies f equal 
to or greater than 0. 
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APPENDIX C 


DEVELOPMENT OF LATTICE MODEL EQUATION 

A ray is traced through the lattice model as depicted in figure 4. 
Starting with the "source" Eq, transmissions and first-order multiple reflec- 
tions at each interface encountered are summed to determine and B2. 

The lattice front surface has the reflection coefficient R as a result 
of coupling to a buffered transducer as in figure 1. There are n interfaces 
within the lattice, each with the identical reflection coefficient G. The 
back surface interfaced with a vacuum (or air) has a reflection coefficient of 
unity. 

In this appendix, as through the main text, the quantities A, B, G, E, 

F, and R are frequency- domain quantities (i.e., Fourier transforms of cor- 
responding time-domain quantities). Reflection coefficients are based on the 
power definition given in appendix A and therefore B^ and B£ must be taken 
as power (energy) spectra of corresponding waveforms. 

Ignoring backscatter (second- order multiple reflections) caused by Eq, 
the first back echo is 


B 1 = V H ~ R) (cl) 

where H is the Heaviside unit step function H = |H(f)| = 1 (appendix B) . 

The second back echo is 



£ G q ~V q 


_1 


2(n-l> 


(C2) 


where A is unspecified intrinsic and/or extrinsic attenuation within the 
lattice layer and is identical in each layer. The echo B2 is obtained by 
summing the internal transmissions and reflections Fi, F2, .... F n and Ej, 
e 2* •••» E n . For forward-going terns (e.g. , F^ to F n ) 
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where t 

F = E q R(H - G) A 

F 12 = E o r2<3(H " G)a3 
f 13 = e q r 3 g 2 (h - G)A 5 

F. = E n R q G q_1 (H - G)A 2q_1 
lq 0 

and 

F = F (H - G) Y G 2(q_1) A 2q_1 

n (n-1) 

For returning echo terms (e.g., E, to E ) 

1 n 

E 1 = F n (H - G) £ G q-1 A 2q = E n + E 12 + . . . + E lq 
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Combining equations (Cl) and (C2) yields 


RB, 


= (H 


G)2 n A 2(n+l) 
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£ G q 'V (q ‘ 1) 
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(CIO) 


In equation (CIO) indicates that the summations are to be multiplied 

together and that additional factors comprised of summations of second-order 
multiple reflections are omitted. These additional factors arise from back- 
transmissions and re-reflections that originate with first-order terms (e.g., 
E 1 , E 2 , .... E n and F lt F 2 , .... F n ) . 

Ignoring second-order multiple reflections within the lattice layers, the 
lattice transfer function T is obtained by combining equation (CIO) with 
equation (3) to obtain equation (16) in the main text. 
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Figure 1. - Schematic of signal acquisition and data-processing stages for determining frequency 
dependence of attenuation coefficient by using broadband ultrasonic pulse-echo method. 
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Fiqure 2. - Diagram of broadband ultrasonic echo system for determining fre- 
quency dependence of attenuation coefficient and deriving material micro- 
structure transfer function. 




Figure 3. - Delta function and Heaviside step function representing time- and frequency-domain coun- 
terparts, respectively, of reflection coefficient, where I Rtf) 1 ■ l Htf) I ■ 1 at free surface (appendix B). 
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Figure 4. - Diagram of lattice model for synthesizing material structure transfer 
function (appendix C). 
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(a) Grain-delay probability distribution function. 
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Frequency, f, MHz 

(b) Magnitude and phase spectra of grain-delay distribution function. 

Figure 5. - Time-domain probability distribution function for grain-delay 
(transit time) and its Fourier transform, magnitude, and phase. 
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(a) Original experimental data and curves for iron. 



0 4 8 12 16 20 

(b) Simulated attenuation data curves for iron (G ■ K - 0.0134, 
h- 0. 00002, v- 0.59 cm /ms. 




Frequency, f, MHz 

(c) Simulated curves for iron to 100 MHz (G - K - 0. 0067, h ■ 0. 00002, (d) Simulated curves for iron to 100 MHz (G * K * 0. 0134, h - 0. 00002, 

v- 0. 59 cm/[is). v- 0, 59 cm /ps). 

Figure 6. - Comparison of experimental attenuation -coefficient-versus-frequency curves for iron with curves generated with equation (34). Oriqinal 
experimental data from Merkulov (1957). 





Attenuation coefficient a, nepers/cm 




(b) Simulated attenuation data curves for copper (G - K ■ 0.0124, 
h- 0.0003, V- 0.47 cm /|is). 


Figure 7. - Comparison of experimental attenuation-coefficient-versus- 
frequency curves for copper with curves generated with equation (34). 
Original experimental data and curves from Generazio (1985). 
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(b) Simulated attenuation data curves for nickel G ■ K * 0. 0144, 
h a 0. 0003, v * 0. 563 cm /ps). 

Figure 8.- Comparison of experimental attenuation-coefficient-versus- 
frequency curves for nickel with curves generated with equation (34). 
Original experimental data and curves from Generazio (1985). 
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(a) Mean grain size, 2 pm. c - 1. 016xl0“ 4 ; m - 1. 835 ; least squares coefficient - 0. 99954. 



log frequency, MHz 


(b) Mean grain size, 4 pm. c - 1.143xl0“ 4 ; m * 1. 975; least squares coefficient - 0. 99947. 

Figure 9. - Simulation of attenuation-coefficient-versus-frequency data in 10- to 100-MHz range using equation (34) and 
determination of c and m values by least squares fit to equation (14) from 15 to 100 MHz. v - o 56 cm /us- K - 0 007* 
G ■ 0. 007 ; h » 2. 0x10 . ‘ * 










(g) Mean grain size, 14 pm. c ■ 1. 487xl0~ 4 ; m - 2. 297; least squares coefficient ■ 0. 99958. 



log frequency, MHz 

(h) Mean grain size, 16 pm. c ■ 1. 729x10"^; m - 2. 301; least squares coefficient - 0. 99971. 

Figure 9. - Continued. 
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(b) 250-grade ma raging steel. 


Figure 10.- Interrelations between experimentally determined 
attenuation parameters c and m from equation (14) for 
two ma raging steels. Data from Vary (1978). 
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Exponent, m 


(c) Parametric field for v - 0,55 cm/ps, a - 10 pm, G c 0. 007. 

Figure 11. - Parametric fields of c versus m based on equations 
(14) and (34) for assigned values of mean grain size a, velocity 
v, hysteresis constant h, average reflection coefficient G, and 
elastic anisotropy factor K. Using the curve-fitting method 
illustrated in figure 9 the fields contain only points for which 
the correlation coefficient exceeds 0. 999 over the frequency 
range 15 to 100 MHz. Data for 250-grade maraging steel from 
figure 10 are coplotted within fields. 
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